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SUMMARY 

The shock standoff distance ahead of a general rounded stagnation 
point has keen estimated under the assumption of a constant-density 
shock layer. It is found that, with the exception of almost-two- 
dimensional bodies with very strong shock waves, the present theoretical 
calculations and the experimental data of Zakkay and Visich for toroids 
are well represented by the relation 

^3D _ / ^ax syrrA / 2 \ 

r s ,x = v r s A k + v 

where A is the shock standoff distance, Rs,x i- s the smaller principal 
shock radius, and K is the ratio of the smaller to the larger of the 
principal shock radii. 


INTRODUCTION 

In recent years much attention has been given to the problem of the 
inviscid flow about blunt shapes, particularly about bodies of revolu- 
tion and cylinders. The theories proposed range from those that con- 

the shock layer to be of constant density to exact numerical inte- 
grations of the compressible flow equations. (A summary of these tech- 
niques and an extensive list of references are given in ref. 1.) The 
inviscid flow about general blunt shapes (finite bodies with unequal 
principal curvatures) has, however, received very little attention. 

Hayes (ref. 2) has derived an expression for the shock standoff distance 
ahead of a general stagnation point for a constant-density shock layer. 
This expression includes centrifugal corrections to the pressure and 
velocity distributions. It was, however, not evaluated in reference 2, 
perhaps because the centrifugal effects depended on shock shape away 
from the stagnation point and a representative three-dimensional shock 
or body shape is difficult to choose. 
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The present analysis retreats from that of Hayes (ref. 2) in that 
centrifugal effects are neglected. In this respect it is the three- 
dimensional analog of an earlier analysis by Hayes (ref. 3). The shock 
location is a function only of the density retio across the shock and 
the ratio of the principal radii of curvature of the shock. It is inde- 
pendent of whether the body is, for example, an ellipsoid or a toroid. 

Because of the many assumptions and simplifications that are made 
in the present analysis, it is not reasonable to expect the theory to 
yield precise absolute values for the shock standoff distance ahead of a 
general stagnation point. It will be shown, however, to give reasonably 
good results for the ratio of this standoff distance to that for an 
axially symmetric body. Then, this ratio together with exact solutions 
for axially symmetric bodies such as those of Van Dyke and Gordon (ref. 
4) may give reasonably good results for general bodies. 

The main analysis is preceded by a very approximate calculation 
(OVERSIMPLIFIED METHOD) which, in spite of its crudity, yields a result 
that is in general agreement with the main analysis. The present re- 
sults are also compared with the experimental data of Zakkay and Visich 
(ref. 5) . 


SYMBOLS 

C shock-layer velocity gradient, eq. (l) 

C x , C z shock-layer velocity gradients in principal directions 

K ratio of smaller to larger principal radius of shock wave. 



k ratio of free- stream density to that behind a normal shock. 

Pco/P, 

K, free- stream Mach number 

p static pressure 

R radius of curvature 

R X ,R Z principal radii of curvature 

u, v,w velocities in x, y, and z directions, respectively 
S/x 


X 
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coordinate in shock layer identified with smaller principal 
radius of curvature of shock 

direction normal to hody 

coordinate in shock layer identified with larger principal ra- 
dius of curvature of shock 

V 2k/ (1 - 2k) 


X ratio of specific heats 

A shock standoff distance 

£ dummy variable in shock surface in z 

| dummy variable in shock surface in x 

p density 

Subscripts: 

ax sym axially symmetric 
b body 

s shock 

1 conditions in shock layer 

2D two - dimens i onal 

3D three-dimensional 

co free-stream conditions 


direction 

direction 


OVERSIMPLIFIED METHOD 


In this approximation the constant-density shock layer is of uni- 
form thickness in the neighborhood of the stagnation point, which is the 
region under consideration. The free-stream Mach number is very large, 

the shock layer is very thin f— « 1, — « 1), and the pressure distri- 

\ R s R s / 

bution is assumed to be Newtonian and constant along a normal to the 
body. The velocity in the shock layer is also assumed constant along a 
normal to the body and equal to the velocity adjacent to the body. 



4 


Within these assumptions the shock-layer velocity may be expressed 
U 1 ~ Cx. R° r ^ 1 an( ^ ^ °°> the velocity gradient is 

[y I T Poo Y* - 1 1 I > 

c = r; But in this limit k = ^ ~ Wr or r - r^-p s 

the shock-layer velocity gradient can be wr: tten 


Uco f 

■ Rs V; 


2k 
1 + k 


Three cases are now considered: 
Case I: Axially Symmetric Flow 


Shock 



(i) 


From continuity considerations (see sketch), 

P„ u m 30 c 2 = P-.U, 2itx A 
oo oo ^11 ax sym 

Using equation (l), the shock-layer thickness becomes 

a ~ JTTT 

^ax sym ~ 

Case II: Two-Dimensional Flow 


The continuity equation per unit span is 

Pco^x = PpUp A 2d 
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so that the shock-layer thickness becomes 

A , ^ /I +" k 
A 2D ~ y 2lT 

which is just twice that for the axially symmetric case. 
Case III: Three-Dimensional Flow 



Let Rs,x he 'the smaller of the principal shock radii and let 
K = R_ Wr o For this case, the continuity equation in each quadrant 

S y S y Z 

can be written 


Poo u co xz = Pl u l z A 3D + P 1 W 1 X A 3D 


or 


P u 

OO 00 


A ® ' n (c x + c z> 




1 + k 
2k 


1 +. 1 


hs , X ®s , z 


In terms of Rs,x and K, 


kR, 




1 + k 
2k 


a 3D = 


1 + K 
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Note that in all cases the standoff distances are of the form 
(VRs,x) ~ which is not appropriate to rounded bodies. This incor- 

rect form results from assuming that at a gj ven x and z the veloc- 
ities in the interior of the shock layer are equal to those at the body 
surface, so that in effect all the fluid in the shock layer entered 
through the normal portion of the shock wave as did that fluid at the 
body surface. Nevertheless, the ratio of the preceding results for a 
given density ratio k, namely. 



will turn out to well represent the results of the following analysis. 

ANALYSIS 

The analysis for the three-dimensional case will now be somewhat 
improved. In particular, the variation of velocities u^ and v^ 
across the shock layer will be taken into account. 

Consider the flow in a quadrant of the constant-density shock layer 
ahead of a three-dimensional stagnation point. The coordinate system is 
considered locally Cartesian. 



n / r 7T 



7 


co 
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00 
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The pressure distribution is taken to be Newtonian. It is o 
function of x and z orF v , lives independent of ^ho n i. 
y. in the neighborhood cl* the stagnation point it may be wr 


■nsldered a 
coordinate 
listen for 


H,c » k « 1 



With the aid of Bernoulli's relation, the velocities at a pcim in the 
shock layer can be estimated. Consider the streamline that passes 
through the shock at the point (?,£)• Bernoulli's relation c eases: 


Pi + P, 


U 1 + V 1 ' W 1 


ll’l + Dp 


u 2 + v 2 + w 2 


f . r 


From the oblique shock relations in the limit M*, » 1, k « 1 


4) 


u i. + v i + W 1 1 , > 


1 - (1 - k 2 ) 1 - 


p2 p 2 


4 * • • 


f M 


Upon substituting relations (3) and (5) into equation (4), 


u 2 + V 2 + w 2 


L 


R. 


s,x 


R. 


s, z 


+ 2k 


fx 2 - * 2 + Z 2 - t; 2 


+ /?(k 2 ) 


R 


s, x 


R5 


O') 


For the portion of the shock layer being considered, the principal 
shock radii are assumed constant. From equation (3), the pressure gra- 
dient in a given principal direction is then independent of the other 
principal coordinate so that Up = up(x,|) and Vp = Wp(z,p). Since, 

in addition, the normal velocity v/u*, is of order k (this can be ver- 
ified a posteriori), the shock-layer velocities in the x and z direc 
tions are, respectively. 





u r 


r> 




+ €( k 2 ) 


W 

U 


2 

1 

2 



+ 0(k 2 ) 


( 7 ) 

(8) 


, U 1 X 

For the stagnation streamline ( £• = £ = o), — « /2k and 

U oo R S,X 

W 1 Z j 

— ^ p — -/2k. These are the velocities obtained upon assuming as in 

U 00 ♦ Z 

the af ©represented OVERSIMPLIFIED METHOD that all the flow enters the 
shock layer through the normal portion of the how shock wave. Currently, 
however, the shock-layer velocity up at depends also on the height 
' ut which the streamline in question traversed the curved how shock 
wave. The velocity wq at z depends similarly on t • 

The shock displacement distance is now to he determined. From con- 
tinuity, the mass flow into the shock-layer quadrant through face ABCD 
must equal the sum of the flows leaving through the top face BCEF and the 
si.de face ABFG. 



But since the normal velocities through the top and side faces depend on 
I and t, respectively, the division of outflow between the top and side 
faces must he determined to properly carry cut the mass balance. In other 
words, the curve BIID must he found such that the inflow through area BCDH 
exits the top face of the control volume while the inflow through ABHD 
exits the side face of the control volume. The line BHD is the inter- 
section of the stream surface through the dividing line BF with the shock 
wave and is found by tracing streamlines hack from line BF, which has the 
coordinates (x,z). Once this is done, the shock displacement distance 
can he ’written either as 



A .Cl O 
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h 




Upon substitution of the velocities from equations (7) and (8) and 
integrating between the prescribed limits 



the equation of line BHD in the form suitable for equation (9b) is found 
to be 



(id 


( 12 ) 


From equations (9b) and (ll) the expression for standoff distance 
is 



( 13 ) 
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By letting X = l/x., equation (13) can be written 


a 3D 



a dX 


sinh 


[(1 - K) 


sinh^a ^ K sinh' 


-1 




X 2 + a 2 


(14) 


Equation (14) is exactly the expression that would be obtained by neg- 
lecting centrifugal effects in equation (4.5.7) of Hayes 1 analysis 
(ref. 2). 

In general, equation (14) must be integrated numerically. However, 
in the special cases of two-dimensional and axially symmetric flow, 
closed-form expressions can be obtained. These are: 

K = 0, two-dimensional flow: 


a 2D = 


kR, 


S, X 


/ 


dX 


I yj x 2 7 a 2 VT^ 2k 


M S,)I . ^/1 - 2k 

-7?=^ sinh y— — 


(15) 


K = 1, axially symmetric flow: 


^ax 


kR, 


S;X 


sym 


/ 


X dX 


kR, 


s,x 


VI - 2k J yjx 2 + a 2 1 + V2k 


;ie) 


These two special results were obtained by Hayes (ref. 3). 


Equations (15) and (16) are now comparer, with exact solutions. In 
figure l(a), it is seen that equation (16) agrees rather closely with 
the calculations of Van Dyke and Gordon (ref. 4) for k <j 0.1, while at 
higher density ratios (k <, 0.3) it overestimates the exact solutions by 
no more than 10 percent. Equation (15) for 1 wo-dimensional flow is less 
satisfactory in that it overestimates the results of Van Dyke (ref. 6), 
Belotserkovskii (ref. 7), and Uchida and Yasihara (ref. 8) by 20 to 50 
percent. This poor agreement is reflected also in the ratio of two- 
dimensional to axially symmetric standoff distances as shown in figure 
i(b). 
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For the three-dimensional case (0 < K < 1) with density ratio k 
greater than zero, equation (14) was numerically integrated on a desk 
calculator using Simpson’s rule. The results are given in table I and 
are also plotted in figure 2. The results shown for k = 0 (table l(b), 
fig. 2(b)) were obtained by evaluating the integrand of equation (14) 
for a -> 0 and then integrating. The result is 




(17) 


The curves of figure 2 are seen to be quite regular, which is per- 
haps not very surprising. However, the resulting ratios of three- 
dimensional to axially symmetric standoff distances given in table II 
and plotted in figure 3 show a more interesting result, namely, that 
most calculated points show surprisingly good agreement with the crudely 
derived expression (2) from the OVERSIMPLIFIED METHOD. The exceptions 
are for bodies approaching the two-dimensional, for example, for 
K < 0.2 at density ratios k < 0-1. 


Before proceeding to a comparison with experiment, it must be real- 
ized that the present analysis yields no information regarding the vari- 
ation of shock-layer thickness about the body and leaves the body shape 
unspecified. It is therefore not suited to determining the shock stand- 
off distance ahead of a given body. Considering some of the results for 
a sphere (ref. 2), it is doubtful whether much is gained by considering 
centrifugal effects and higher order pressure terms in a constant-density 
approach. It seems rather that, if the constant-density solution is to 
be improved on, the compressible flow equations should be solved exactly 
as done in references 4, 6, 7, and 8 for axially symmetric and two- 
dimensional bodies. 


COMPARISON WITH EXPERIMENT 

The only pertinent experiments known to the author are those at 
Moo = 3 and M» = 8 by Zakkay and Visich (ref. 5). The three-dimensional 
body tested was a toroid. The vital statistics of the experiments and the 
theoretical comparison are given in the following table (shock radii un- 
fortunately had to be measured from the schlieren photographs presented 
in ref. 5) . 
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H, = 3 

Hx, = 8 

Experiment - (ref. 5) 



^b,x/R*b, z 

0 255 

0.255 

^ = ^s,x/^s,z 

■ 45 

.34 

A / R b^x 

408 

.267 

A/^s,X 

194 

.152 

Theory 



k 

0-259 

0.180 

4/R s ,x ( from 2 ) 

. 208 

.169 

(eg. 18) 

187 

.157 

(eg. 19) 

■ 186 

.156 


The experimental data for A/R s ^ x are compared with three theoretical 

estimates. The first is that value taken directly from figure 2(a) for 
the pertinent density ratio and shock radius ratio. The second calcula- 
tion is according to the relation 



" ( A 3D \ 

^3D ( ^ax sym\ 

\ R s,x/ 

Rs ’ x ~ \ *= /van Efyke 

A^ax sym\ 

Lv R s ) 


where the ratio of three-dimensional to axially symmetric standoff dis- 
tances for the proper density ratio k is taken from figure 3. The 
third calculation uses equation (2) for the aforementioned ratio: 


A5D 
Rs, x 




(19) 


For both equations (18) and (19) the axially symmetric standoff distance 
is taken from Van Dyke and Gordon (ref. 4) * whose data are partially 
shown in figure 1. 


The theoretical estimates agree well wiuh the experimental standoff 
distances; in fact, the agreement is better ahan might be expected, con- 
sidering the author's ability to measure shock radii from photographs. 
The results using equation (18) or (19) are ^dthin 5 percent of the 
measured values. The theoretical estimate from figure 2 is somewhat 
higher, which indicates primarily that the present constant-density ap- 
proximation overestimates the absolute standoff distance for all bodies. 
This has already been shown for two-dimensional and axially symmetric 
bodies in figure 1. 
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CONCLUDING REMARKS 

The shock standoff distance ahead of a general rounded stagnation 
point has been estimated under the assumption of a constant-density 
shock layer. It is found that many of the present theoretical calcula- 
tions as well as the two experimental points of Zakkay and Visich are 
well represented by the relation 



A where K is the ratio of the smaller to the larger principal shock ra- 

dius. The exceptional cases are bodies approaching the two-dimensional 
(K <, 0.2) with shock layers whose density is much larger than that of 
the free stream (k < 0.10). In comparing with experiment, the axially 
symmetric standoff distance was taken from the exact solutions of Van 
Dyke and Gordon. 

Unfortunately, a constant-density theory gives no information re- 
garding the body shape corresponding to a given shock wave and is there- 
fore not suited for obtaining the standoff distance ahead of a given 
body. In fact, from the experience with the two-dimensional and axially 
symmetric problems, about the only solutions that adequately relate the 
flow field to the body are the exact compressible flow solutions. The 
present results may nevertheless serve as a semiquantitative guide to 
the phenomenon in the absence of exact solutions. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, June 5, 1961 
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standoff distances. 

Figure 1. - Comparison of exact solutions for two- 
dimensional and axially symmetric shock standoff 
distances with Hayes' simple theory (ref. 3). 
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( a ) a 3d/ r ;;; ,x- 


Figure 2. - Shock standoff distance for three-dimensional bodies. 






Figure 3. - Ratio of three-dimensional to axially symmetric 
shock standoff distances. 


NASA - Langley Field, Va. E-1278 




